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Abstract 

The perovskite-like compounds AC113R114O12 (A=Na,Ca,Sr,La,Nd) are studied by 
means of density functional theory based electronic structure calculations using the 
augmented spherical wave (ASW) method. The electronic properties are strongly 
influenced by covalent type bonding between transition metal d and oxygen p states. 
The characteristic tilting of the RuC"6 octahedra arises mainly from the Cu-0 bond- 
ing, allowing for optimal bond lengths between these two atoms. Our results provide 
a deeper understanding of octahedral tilting as a universal mechanism, applicable 
to a large variety of multinary compounds. 
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1 Introduction 



Despite its striking simplicity the perovskite structure, ABX3, contains nu- 
merous crystallographic variations giving room for a huge class of compounds 
[1,2]. Among these the vast majority comprises oxides and fluorides [3], but 
also chlorides, hydrides, oxynitrides, and sulfides. The great interest in these 
materials is motivated by a variety of exciting dielectric, magnetic, electrical, 
optical, and catalytic properties, with several technological applications. 
For optimal tailoring of materials much work on the perovskite-related com- 
pounds concentrates on the interrelations between deviations from the ideal 
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perovskite structure and physical properties. These crystallographic devia- 
tions can be grouped into three different mechanisms [1,2]: While distortions 
of the characteristic BX 6 octahedra and cation displacements within the oc- 
tahedra are mainly driven by electronic instabilities of the octahedral metal 
ion, octahedral tiltings represent the most common deviation and are observed 
when the A cation is too small for filling the space between regularly ordered 
octahedra. In this situation, octahedral tilting represents the lowest energy 
distortion, since it allows to adjust the A-0 distances while leaving the first 
coordination sphere of the M cation intact, and changing mainly the soft M- 
O-M bond angle rather than the M-0 bond length. 

Of particular interest in this context are the perovskite-like compounds of 
the kind AA3B4O12, where the BOB bond angle is the relevant quantity 
in order to optimize both the A-0 and A'-O bond lengths. Studying these 
materials thus promises deeper insight into the mechanisms driving the oc- 
tahedral tilting. In addition, this material class shows extraordinary physical 
properties. While CaCu3Mn40i2 is a ferromagnetic semiconductor with high 
Curie temperature and large magnetoresistance [4], CaCusTi40i2 is known for 
its unusually high low-frequency dielectric constants [5,6,7]. In contrast, the 
metallic ruthenates of the type ACu 3 Ru 4 0i2 (A=Na,Ca,La) are Pauli para- 
magnets [8] and display valence degeneracy [7]. The crystal structures of the 
latter class have been studied recently by EXAFS and XRD, revealing consid- 
erable deviations of the measured interatomic distances from those expected 
from a bond valence approach [9]. 

The present study investigates the relationship between the octahedral tiltings 
and the electronic properties of the ruthenates ACu 3 Ru 4 0i2 (A=Na,Ca,Sr,La,Nd). 
The discussion is based on electronic structure calculations within density 
functional theory, which are known for their predictive power (see e.g. Refs. 
[10,11] and references therein). As a result, the specific bond lengths as well 
as the discrepancies between the experimental interatomic distances of the 
class ACu 3 Ru 4 0i2 and those proposed by the bond valance approach can be 
understood from the electronic properties and the type of chemical bonding. 



2 Crystal Structure 

The compounds of the class AA3B4O12 crystallize in a body-centered cubic 
2x2x2 superstructure of the simple cubic perovskite structure with space 
group Im3 (No. 204). The atoms are located at the following Wyckoff posi- 
tions: A (2a) at (0,0,0), A' (66) at (|, 0, 0), B (8c) at (±, \, \), and O (24#) at 
(x,y,0). The crystal structure is shown in Fig. 1, where the BOg octahedra 
are easily identified. The tilting of the octahedra is controlled by the oxygen 
position, which, for all the ruthenates under study, is specified by x Q ps 0.175 
and yo ps 0.305. The deviation from the ideal values xq = yo = 0.25 causes a 
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Fig. 1. Crystal structure of AA3B4O12. Large and medium size spheres denote A- 
and Cu-atoms, respectively, while octahedra represent the BOg units. 

rotation of the octahedra around the [111] axis. As a consequence, the A sites 
are partitioned into the A site at (0,0,0) and the A' site, which takes three 
quarters of the original A positions. In particular, the A' sites are strongly 
affected by the oxygen atom shifts away from the ideal position and thus their 
surrounding experiences the largest distortion. As a consequence, three dif- 
ferent A'-O distances of « 2.0 A, 2.8 A, and 3.2 A appear, each forming an 
approximately square-planar coordination and thus providing an ideal geom- 
etry for a Jahn- Teller active ion as e.g. Cu 2+ [9]. The Cu0 4 squares can be 
easily identified in Fig. 1. 

For the ruthenates ACU3RU4O12 (A=Na,Ca,Sr,La,Nd) studied in the present 
work, crystallographic data were determined by Ebbinghaus et al. [9]. The 
lattice constants are similar to those reported by Labeau et al. [8] as well as 
Subramanian and Sleight [7]. Both the lattice constants and oxygen parame- 
ters of these compounds show only slight deviations and amount to a ph 7.4 A 
as well as x Q ~ 0.175 and yo ~ 0.305 as mentioned above. 



3 Method of Calculation 



The electronic structure calculations were performed within the framework of 
density functional theory and the local density approximation using the aug- 
mented spherical wave (ASW) method [12,13]. In order to represent the cor- 
rect shape of the crystal potential in the large voids, additional augmentation 
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spheres were inserted into the open crystal structures. Optimal augmentation 
sphere positions as well as radii of all spheres were found automatically by the 
sphere geometry optimization (SGO) algorithm [14]. The basis sets comprised 
Cu 4s, Ap, 3d, Ru 5s, 5p, Ad, Af, and O 2s, 2p as well as states of the additional 
augmentation spheres. In addition, Na 3s, 3p, 3d; Ca 4s, Ap, 3d, Af; Sr 5s, 
5p, Ad, Af; La/Nd 6s, Qp, 5d, Af states were used for the A cation. The Bril- 
louin zone sampling was done using an increasing number of k points within 
the irreducible wedge ranging from 17 to 1255 points to ensure convergence 
of the results with respect to the k-space grid. In addition to the analysis of 
the (partial) densities of states we will also address the question of chemical 
bonding in terms of the covalent bond energy, which was recently proposed 
[15] as an extension to the concept of the crystal orbital overlap population 
(COOP) [16]. In short, the covalent bond energy covers those contributions 
to the total energy of a crystal, which arise from the hybridization of orbitals 
located at different atoms. 



4 Results and Discussion 

Partial densities of states (DOS) for all five compounds under consideration 
are displayed in Fig. 2. The electronic structure in the energy window shown is 
dominated by Cu 3d, Ru Ad, and O 2p states. In addition, Nd Af states cause 
a sharp peak at and slightly above E F for NdCu 3 Ru 4 0i2. The gross features 
are very similar for all compounds. Three energy ranges can be distinguished. 
In the broad interval from about -8 to -3eV oxygen 2p states dominate but 
are complemented by considerable contributions from the transition metal d 
states due to covalent bonding. Whereas in the energy range from -3 to -2 eV a 
sharp Cu 3d peak appears, electronic states between -2 and 0.5 eV, and hence 
metallic conductivity, derive mainly from broader Ru Ad bands. Covalent type 
bonding leads to finite oxygen 2p contributions above -3 eV. This is confirmed 
by covalent bond energy curves calculated for the sodium compound as shown 
in Fig. 3. Both the Cu-0 and the Ru O curves are negative and positive at 
energies below and above -3 eV, respectively, indicative of bonding and anti- 
bonding states in the energy regions, where the O 2p and the transition metal 
d states dominate. However, due to the higher d electron count of Cu as com- 
pared to Ru, the Cu-0 antibonding states are filled, leading to a vanishing 
contribution of these bonds to the chemical stability. In contrast, the inte- 
grated Ru-0 covalent bond energy curve at E F has a finite value indicative 
of the stabilizing net contribution of this bond. Finally, metal-metal bonding 
results in small peaks between -3 and -2 eV and above the Fermi energy. 
All states which are not displayed in Fig. 2 play only a very small role in the 
energy interval shown. In particular, apart from the Nd Af states, the A cation 
electrons do not contribute to the covalent-type bonding but are distributed 
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Fig. 2. Partial Cu 3d, Ru 4d, and O 2p densities of states (DOS) of ACU3RU4O12 
per formula unit. 4/ states have been included for the Nd-compound. 



over the solid. As a consequence, changes at the A cation site show up mainly 
in the modified total electron count. 

Ebbinghaus et al. [9] related the measured atomic distances for ACU3RU4O12 
(A=Na,Ca,Sr,La,Nd) to the predictions of the bond valence model. The latter 
is an empirical model, which connects the valences and the mutal distances 
of the atoms [17]. It starts out from ionic configurations but looses accuracy 
in covalently bonded crystals. For the divalent A cations the Ru-0 distances 
were found to agree with the expected values. In the case of A=Na the bond 
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Fig. 3. Covalent bond energy curves for NaCu3Ru40i2- The corresponding inte- 
grated covalent bond energies are given in the inset. 



length is longer, while for A=La,Nd it is shorter than predicted by the bond 
valence model. 

In order to understand the influence of octahedral tilting on the electronic 
properties and to analyze the mentioned discrepancies reported by Ebbing- 
haus et al. we performed additional calculations using hypothetical crystal 
structures. In these setups the oxygen atoms were shifted away from the mea- 
sured positions to the ideal position x — 0.25, yo = 0.25 and to x Q = 0.15, 
yo = 0.35. These positions correspond to zero tilting of the octahedra and 
to an exaggerated tilting, respectively. Note that shifting the oxygen position 
mainly affects the Cu-0 bond lengths, while the effect for the Ru-0 and A-0 
distances is of second order. The partial DOS for the hypothetical structures 
of NaCu3Ru40i2 are displayed in Fig. 4. Considerable changes as compared 
to the results for the experimental crystal structure (see Fig. 2) are observed. 
Turning to the ideal structure first, we observe (i) a considerable broadening 
of the O 2p dominated states and (ii) a striking sharpening of the Cu 3d states. 
Both effects are readily understood from the oxygen shifts. The latter cause 
an increase of the Cu-0 distances and, to a lesser degree, a decrease of the 
Ru-0 distances, leading to a slightly larger 2p-Ad overlap. As a consequence, 
larger Ru Ad contributions at energies below -3 eV are observed, whereas the 
Cu 3d fraction is almost negligible in this region. At the same time, Ru Ad 
and O 2p admixture to the sharp Cu 3d peak is likewise very small. To sum 
up, shifting the oxygen atoms to their ideal positions leads to an effective 
decoupling of the compound into Ru0 6 octahedra and Cu atoms, which are 
effectively separated and thus do no longer form a stable solid. 
On the contrary, for the hypothetical structure with the oxygen shift away 
from the ideal position exaggerated (i) the oxygen states have moved to lower 
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Fig. 4. Partial Cu 3d, Ru 4d, and O 2p densities of states (DOS) of hypothetically 
distorted variants of NaCusRu40i2 per formula unit. 

energies and (ii) the Cu 3d states do no longer form a distinct peak but have 
attained a rather itinerant character. The latter goes along with an increased 
similarity of the Cu 3d and O 2p partial DOS. The increased covalent type 
bonding between these two states is also reflected by large contributions to the 
covalent bond energy (not shown), negative below -4eV and positive above, 
leading to a net antibonding integral at Ep. The increased 2p-3d bonding is 
in contrast to a slight decrease of the bonding between the Ru 4d and O 2p 
states due to the larger extension of the Ru0 6 octahedra. Thus the hypothet- 
ical structure with exaggerated oxygen shifts likewise suffers from an increase 
of antibonding at the expense of bonding states. 

Together with the instability of the hypothetical structure with ideal oxygen 
positions we conclude that the experimental structure is due to an optimal 
balance of different types of bonds, namely 2p-3d and 2p-4d. This balance is 
less affected by the A cations, which do not take part in covalent bonding but 
have lost their outer electrons. For this reason, from the point of view of chem- 
ical bonding, there is no means to optimize the A-0 bond length. Moreover, 
as for the Ru O distance, the A-0 distance is only to a small degree affected 
by the tilting. Finally, the different types of bonding present in these com- 
pounds give an additional explanation for the deviations of the experimental 
interatomic distances from those expected from the bond valence model. 



5 Conclusions 

The electronic properties of the perovskite-like ruthenates ACu 3 Ru 4 0i2 are 
governed by strong covalent bonding between the transition metal d and the 
oxygen 2p electrons. While Ru-0 bonds via the corresponding bond-lengths 
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influence mainly the size of the octahedra and hence the lattice constant, the 
octahedral tilting is to a large part due to the Cu-0 bonding. 
As compared to the ideal perovskite structure the tilting of the characteristic 
Ru0 6 octahedra has only little effect on the Ru-0 bonding. In contrast, the 
Cu-0 distances undergo serious changes, leading eventually to the square- 
planar coordination of copper. Finally, the unusual A-0 distances and the 
failure of the bond valence approach are due to the fact that their response 
to the tilting is only small and the structure does not offer any possibility to 
simultaneously optimize the A-0 and the Cu-0 bond lengths. 
In conclusion, octahedral tilting as discussed for ACu 3 Ru 4 0i2, arises from 
the interplay of covalent bonding between different atomic species, leading to 
optimized bond lengths between all atoms involved. This general principle is 
of rather universal nature and applies to a large variety of compounds. 
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